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Abstract 
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1. FACTORIZATION &: THE NLO MODEL 

A generic inclusive cross section for the process A + 
B ^ F + X with observed final-state system F, of 
total mass Q, can be expressed as 

4 daAB^FX , , 2\ ^ J / 2\ 
^ (Pa/A{Xa,fJ- ) ® (pi, / B^Xb , ) 

® ^ab^FX {z,Q,fl) , (1) 

with z = Q^/xaXbS. The aab are partonic hard- 
scattering functions, a — (TBorn + ioisil^'^)/''^)o'^^'^ + . . . . 
They are known to NLO for most processes in the stan- 
dard model and its popular extensions. Corrections be- 
gin with higher, uncalculated orders in the hard scat- 
tering, which respect the form of Eq. ([T]) . The discus- 
sion is simplified in terms of moments with respect to 

CTAB^FX = / dr T^^^ daAB^Fx/dQ^ 
Jo 

= J2 ^a/A{N,^i^) ^ab~>Fx{N,Q,^i) 4>b/ B {N , , (2) 

where the moments of the (/)'s and o'ab^FX are defined 
similarly. 

Eqs. (|l|) and (||) are starting-points for both the de- 
termination and the application of parton distribution 
functions (PDFs), ifii/n, using 1-loop a's [ |^, |[ ^ 
We may think of this collective enterprise as an "NLO 
model" for the PDFs, and for hadronic hard scattering 
in general. For precision applications we ask how well 
we really know the PDFs [ |, |, ^j. Partly this is a 
question of how well data constrain them, and partly 
it is a question of how well we could know them, given 
finite-order calculations in Eqs. (|l|) and (Q). We will 
not attempt here to assign error estimates to theory. 
We hope, however, to give a sense of how to distin- 
guish ambiguity from uncertainty, and how our partial 
knowledge of higher orders can reduce the latter. 
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2. UNCERTAINTIES, SCHEMES & SCALES 

It is not obvious how to quantify a "theoretical un- 
certainty" , since the idea seems to require us to esti- 
mate corrections that we haven't yet calculated. We 
do not think an unequivocal definition is possible, but 
we can try at least to clarify the concept, by consid- 
ering a hypothetical set of nucleon PDFs determined 
from DIS data alone [ ^ . To make such a determina- 
tion, we would invoke isospin symmetry to reduce the 
set of PDF's to those of the proton, t/fo/p, and then 
measure a set of singlet and nonsinglet structure func- 
tions, which we denote F'-'^K Each factorized structure 
function may be written in moment space as 

(iV, Q)^J2 ^a/p{N, m') , (3) 

a 

in terms of which we may solve for the parton distri- 
butions by inverting the matrix C, 

i 

With "perfect" F's at fixed Q, and with a specific 
approximation for the coefficient functions, we could 
solve for the moment-space distributions numerically, 
without the need of a parameterization. In a world 
of perfect data, but of incompletely known coeffi- 
cient functions, uncertainties in the parton distribu- 
tions would be entirely due to the "theoretical" uncer- 
tainties of the C"s: 

5~^a,p{N,^i)^Y.^^~^^a{N,Q,^x) F^^{N,Q). (5) 

i 

Our question now becomes, how well do we know the 
C's? In fact this is a subtle question, because the coef- 
ficient functions depend on choices of scheme and scale. 

Factorization schemes are procedures for defining co- 
efficient functions perturbatively. For example, choos- 
ing for F2 the LO (quark) coefficient function in Eq. 
(^) defines a DIS scheme (with C independent of /i, 
which is then to be taken as Q in 0). Computing the 
C's from partonic cross sections by minimal subtrac- 
tion to NLO defines an NLO MS scheme, and so on. 
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Once the choices of C"s and fi are made, the PDF's are 
defined uniquely. 

Evolution in an MS or related scheme, enters 
through 

^i-^4>a/HiN,^I^) = -^ab{N,as{^i^))k/HiN,^I^) 

fi-^Ci'HN,Q,fi) = cf{N,Q,^i)Tac{N,as{^Ji'')). (6) 

In principle, by Eq. (^, the scale-dependence of the 
ci*-* exactly cancels that of the PDFs in Eq. and, 
by extension, in Eq. (Q). This cancelation, however, 
requires that each C and the anomalous dimensions T 
be known to all orders in perturbation theory. 

To eliminate ^-dependence up to order we 
need d to order a" and the Tab to a"+^. One- loop 
(NLO) QCD corrections to hard scattering require two- 
loop splitting functions, which are known. The com- 
plete form of the NNLO splitting functions, is still 
somewhere over the horizon [ |^. Even when these 
are known, it will take some time before more than a 
few hadronic hard scattering functions are known at 
NNLO. 

We can clarify the role of higher orders by relating 
structure functions at two scales, Qo ^-nd Q. Once we 
have measured F{N,Qo), we may predict F{N,Q) in 
terms of the relevant anomalous dimensions and coef- 
ficient functions by 



FiN,Q) 



F(iV,Qo) eJ«« ^ 
CiN,Q,Q) 



CiN, Qo,Qo) 



(7) 



This prediction, formally independent of PDFs and in- 
dependent of the factorization scale, has corrections 
from the next, still uncalculated order in the anoma- 
lous dimension and in the ratio of coefficient functions. 
The asymptotic freedom of QCD gives a special role to 
LO: only the one-loop contribution to F diverges with 
Q in the exponent, and contributes to the leading, log- 
arithmic scale breaking. NLO corrections already de- 
crease as the inverse of the logarithm of Q, NNLO 
as two powers of the log. Thus, the theory is self- 
regulating towards high energy, where dependence on 
uncalculated pieces in the coefficients and anomalous 
dimensions becomes less and less important. 

The general successes of the NLO model strongly 
suggest that relations like (Q) are well-satisfied for a 
wide range of observables and values of N (or x) in 
DIS and other processes. This does not mean, however, 
that we have no knowledge of, or use for, information 
from higher orders. In particular, near a; = 1 PDFs are 
rather poorly known [ At the same time, the ratio 



of C's depends on N, and if as IniV is large, it becomes 
important to control higher-order dependence on IniV. 
This is a task usually referred to as resummation, to 
which we now turn. 

3. RESUMMATION 

Let us continue our discussion of DIS, describing 
what is known about the A^-dependence of the coeffi- 
cient functions C, as a step toward understanding the 
role of higher orders. Specializing again for simplic- 
ity to nonsinglet or valence, the resummed coefficient 
function may be written as [ p| 

C'^'^^N, g, ^i) = (7f„LO(A^, Q, + s e^^ois(w,Q,M)^ (g) 

where "sm6" implies a subtraction on C^^'^ to keep 
C'^'^^ exact at order as , and where Cf^^ corresponds to 
the NLO A'-independent ("hard virtual") terms. The 
exponent resums logarithms of A: 



^dis(^,Q,m) 



(9) 
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with A^ = A^e'^^ , and with 
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Eq. (pli| ) is accurate to leading (LL) and next-to-leading 
logarithms (NLL) in A" in the exponent: a™ In' 



m+l 



N 



and a™ In™ A^, respectively. The A^ dependence of the 
ratio C'J°"(Af, Q, Q)/C^^'^In, Q, Q) is shown in Fig. |], 
with Q2 = 1^ 5^ 100 GeV^ At A^ = 1 the ratio is 
unity. It is less than unity for moderate A^, but then 
begins to rise, with a slope that increases strongly for 
small Q. At low and large A^, higher orders can 
be quite important. What does this mean for PDFs? 
We can certainly refit PDFs with resummed coefficient 
functions, and we see that the high moments of such 
PDFs are likely to be quite different from those from 

NLO fits. 

To get a sense of how such an NLL/NLO-MS scheme 
might differ from a classic NLO-MS scheme, we resort 
to a model set of resummed distributions, determined 
as follows. We define valence PDFs in the resummed 
scheme by demanding that their contributions to F2 
match those of the corresponding NLO valence PDFs 
at a fixed Q — Qo, which is ensured by 



Ires/ AT INLO/Ar (-^i *3o : Qo) 



(11) 
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Figure 1. Ratio of Mellin-A^ moments of resummed 
and NLO MS-scheme quark coefficient functions for 
F2 ■ The numbers denote the value of in GeV^ . We 
have chosen fi — Q. 



Using the resummed parton densities from Eq. (pi]), 
we can generate the ratios Q)/F^^'^{x, Q). 

The resuh of this test, picking Qq = 100 GeV^ is 
shown in Fig. |^, for the valence ^2(2;, Q) of the proton, 
with X = 0.55, 0.65, 0.75 and 0.85. The NLO distribu- 
tions were those of [ 0, and the inversion of moments 
was performed as in . The effect of resummation is 
moderate for most Q. At small values of Q, and large x, 
the resummed structure function shows a rather sharp 
upturn. One also finds a gentle decrease toward very 
large Q [ |l^ . We could interpret this difference as the 
uncertainty in the purely NLO valence PDFs implied 
by resummation. 

From this simplified example, we can already see 
that the use of resummed coefficient functions is not 
likely to make drastic differences in global fits to PDFs 
based on DIS data, at least so long as the region of 
small Q^, of 10 GeV^ or below, is avoided at very large 
X. At the same time, it is clear that a resummed fit 
will make some difference at larger x, where PDFs are 
not so well known. We stress that a full global fit will 
be necessary for complete confidence. 

4. RESUMMED HADRONIC SCATTERING 

Processes other than DIS play an important role 
in global fits, and in any case are of paramount phe- 
nomenological interest. Potential sources of large cor- 
rections can be identified quite readily in Eq. (0). At 
higher orders, factors such as as In^iV, can be as large 
as unity over the physically relevant range of z in 
some processes. In this case, they, and their scale de- 
pendence can be competitive with NLO contributions. 
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Figure 2. Ratio of the valence parts of the resummed 
and NLO proton structure function F2(x,Q'^), as a 
function of for various values of Bjorken-x. For 
F2°^, the 'resummed' parton densities have been de- 
termined through Eq. (pT|). 



Since they make up well-defined parts of the correction 
at each higher order, however, it is possible to resum 
them. To better determine PDFs in regions of phase 
space where such corrections are important, we may 
incorporate resummation in the hard-scattering func- 
tions that determine PDFs. 

The Drell-Yan cross section is the benchmark for the 
resummation of logs of 1— z, or equivalently, logarithms 
of the moment variable iV [ ^ , 

+0{l/N). (12) 
The exponent is given in the MS scheme by 

+2^ ^A{as{f^"))\n((^), (13) 

with A as in Eq. (0), and where we have exhibited 
the dependence on the factorization scale, setting the 
renormalization scale to Q. Just as in Eq. ( p^ for 
DIS, Eq. ( p^ ) resums all leading and next-to-leading 
logarithms of N. 

It has been noted in several phenomenological appli- 
cations that threshold resummation, and even fixed- 
order expansions based upon it, significantly reduce 
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sensitivity to the factorization scale [ [T^ . To see why, 
we rewrite the moments of the Drell-Yan cross section 
in resummed form as 



a%{N,Q) 



2^ cj)g/AiN,fi)a^^iN,i 



CTBorn 



DY 



X (j)g/B{N, ^i) e^°^(^>Q'^)/2 + 0{1/N) . (14) 

The exponentials compensate for the InA^ part of 
the evolution of the parton distributions, and the /x- 
dependence of the resummed expression is suppressed 
by a power of the moment variable, 

['^9M(^''^) e^--(~-«-'^)/2 ] =0(1/7V). (15) 

This surprising relation holds because the function 
A(as) in Eq. ( p^ ) equals the residue of the 1/(1 — a;) 
term in the splitting function Pgg. Thus, the remaining 
iV-dependence in a resummed cross section still begins 
at order a^, but the part associated with the 1/(1 — 
term in the splitting functions has been canceled to 
all orders. Of course, the importance of the remain- 
ing sensitivity to fi depends on the kinematics and the 
process. In addition, although resummed cross sec- 
tions can be made independent of /i for all IniV, they 
are still uncertain at next-to-next-to leading logarithm 
in N, simply because we do not know the function A at 
three loops. Notice that none of these results depends 
on using PDFs from a resummed scheme, because MS 
PDFs, whether resummed or NLO, evolve the same 
way. The remaining, uncanceled dependence on the 
scales leaves room for an educated use of scale-setting 
arguments [|l4j. The connection between resummation 
and the elimination of scale dependence has also been 
emphasized in [ |l5| . 

Scale dependence aside, can we in good conscience 
combine resummed hard scattering functions in Eq. (|l|) 
with PDFs from an NLO scheme? This wouldn't make 
much sense if resummation significantly changed the 
coefficient functions with which the PDFs were orig- 
inally fit. As Fig. ^ shows, however, this is unlikely 
to be the case for DIS at moderate x. Thus, it makes 
sense to apply threshold resummation with NLO PDFs 
to processes and regions of phase space where there 
is reason to believe that logs are more important at 
higher orders than for the input data to the NLO fits. 

At the same time, a set of fits that includes thresh- 
old resummation in their hard-scattering functions can 
be made [0, and their comparison to strict NLO fits 
would be quite interesting. Indeed, such a compari- 
son would be a new measure of the influence of higher 



orders. A particularly interesting example might be 
to compare resummed and NLO fits using high-px jet 
data [||. 

5. POWER-SUPPRESSED CORRECTIONS 

In addition to higher orders in as(/x^), Eq. (|^) has 
corrections that fall off as powers of the hard-scattering 
scale Q. In contrast to higher orders, these corrections 
require a generalization of the form of the factorized 
cross section. Often power corrections are parameter- 
ized as h{x)/[{l — x)Q'^] in inclusive DIS, where they 
begin at twist four. In DIS, this higher twist term infiu- 
ences PDFs when included in joint fits with the NLO 
and NNLO models, and vice- versa [ |l|, |l||. As 
in the case with higher orders, such "power-improved" 
fits should be treated as new schemes. 

6. CONCLUSIONS 

The success of NLO fits to DIS and the studies of 
resummation above suggest that over most of the range 
of X, theoretical uncertainties of the NLO model are 
not severe. At the same time, to fit large x with more 
confidence than is now possible may require including 
the resummed coefficient functions. 

Resummation is especially desirable for global fits 
that employ a variety of processes, such as DIS and 
high-pT jet production, which differ in available phase 
space near partonic threshold. In a strictly NLO ap- 
proach, uncalculated large corrections are automati- 
cally incorporated in the PDFs themselves. As a re- 
sult, the NLO model cannot be expected to fit simulta- 
neously the large- a; regions of processes with differing 
logs of 1 — a; in their hard-scattering functions, unless 
these higher-order corrections are taken into account. 

The results illustrated in the figures suggest that 
these considerations may be important in DIS with 
below a few GeV^ and at large x, where they may 
have substantial effects on estimates of higher twist in 
DIS. In hadronic scattering, large-A^ (a; 1) resum- 
mation, which automatically reduces scale dependence, 
may play an even more important role than in DIS. 
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